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Abstract 

Using both numerical and analytical tools we study various fea- 
tures of static, spherically symmetric solutions of the Einstein- Vlasov 
system. In particular, we investigate the possible shapes of their mass- 
energy density and find that they can be multi-peaked, we give numer- 
ical evidence and a partial proof for the conjecture that the Buchdahl 
inequality sup r>0 2m(r)/r < 8/9, m(r) the quasi-local mass, holds for 
all such steady states — both isotropic and anisotropic — , and we give 
numerical evidence and a partial proof for the conjecture that for any 
given microscopic equation of state — both isotropic and anisotropic — 
the resulting one-parameter family of static solutions generates a spiral 
in the radius-mass diagram. 



1 Introduction 

When studying the properties of non- vacuum spacetimes in general relativity 
the choice of the matter model is important. In the present paper matter 
is described as a large ensemble of particles which interact only via the 
gravitational field created by the particles themselves and not via collisions 
between them. The distribution of the particles on phase space is given 
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by a distribution function /. This function satisfies the Vlasov equation, 
a first order conservation law the characteristics of which are the geodesies 
of the spacetime metric, and on the other hand / gives rise to macroscopic 
quantities such as mass-energy density, pressure, and mass current which act 
as source terms in the Einstein field equations. For an introduction to kinetic 
theory in general relativity and the Einstein- Vlasov system in particular we 
refer to p] and [IB] . 

One distinguishing feature of the Einstein- Vlasov system is that it allows 
for a wide variety of static solutions. The purpose of the present paper 
is to analyze spherically symmetric static solutions both numerically and 
analytically. After some general discussion of such steady states in the next 
section and a brief discussion of our numerical approach we investigate three 
different features of the steady states. In Section [3] the energy density as 
a function of the area radius is characterized into three distinct classes. In 
particular, we study steady states which as a function of the area radius are 
supported in an interval [Rq, R\] with Rq > 0, which we call shells, as well as 
states supported in [0, Ri], which we call non-shells. One observation is that 
anisotropic steady states can have a quite rich structure, in particular when 
compared with their Newtonian analogues, the steady states of the Vlasov- 
Poisson system. For instance, the energy density can have arbitrarily many 
peaks and in the case of shells these can be separated by vacuum regions so 
that the steady states consist of rings of Vlasov matter. 

Then we consider the question whether there is an upper bound strictly 
smaller than 1 of the quantity 2M/R\, where M is the ADM mass. For 
isotropic perfect fluid solutions which satisfy the hypothesis that the energy 
density is non-increasing outwards the inequality 

2M 8 , , 

ft < 9' (L1) 

was shown by Buchdahl [7] . We prove a pointwise version of this inequality, 
cf. inequality (14. 5p below, for a class of steady states of the Einstein- Vlasov 
system which includes, but is not restricted to, isotropic states. Moreover, 
we give numerical evidence for the conjecture that all spherically symmetric 
steady states of the Einstein- Vlasov system satisfy this inequality, even if 
the energy density is not radially non-increasing. We mention in passing 
that the inequality is sharp in the case of shells, cf. [2]. 

If one prescribes the microscopic equation of state, i.e., the way in which 
/ depends on the local energy and angular momentum, one obtains a one 
parameter family of steady states. In the last section we investigate the 
relation of the ADM mass to the outer area radius of the steady state along 
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such one-parameter families. We prove the existence of mass-radius spirals 
for the isotropic case by reducing it to a result of Makino [11] for isotropic 
fluids, and we give numerical evidence that these radius-mass spirals are 
present for any microscopic equation of state, independently of isotropy, a 
feature which is again in sharp contrast to the Newtonian situation. 

2 Static solutions of the spherically symmetric 
Einstein-Vlasov system 

In Schwarzschild coordinates the metric of a static spherically symmetric 
spacetime takes the form 

ds 2 = - e Mr) dt 2 + e 2X(r) dr 2 + ^2 + sin 2 

where r > 0, € [0, 7r], cp E [0, 2tt]. Asymptotic flatness is expressed by the 
boundary conditions 

lim A(r) = lim fi(r) = 0, 

r—foo r^oo 

and a regular centre requires 

A(0) = 0. (2.1) 

The static Einstein-Vlasov system is given by the Einstein equations 

e - 2X (2rX' - 1) + 1 = 8vrr 2 p, (2.2) 
e" 2A (2V + 1) - 1 = 8irr 2 p, (2.3) 

together with the static Vlasov equation which can be written as 

w d r f - ((1 + w 2 + L/r V - L/r 3 ) d w f = 0. (2.4) 

The variables w and L can be thought of as the momentum in the radial 
direction and the square of the angular momentum respectively. The matter 
quantities are given by 



/OO /"OO 
/ ^l + w 2 + L/r 2 f(r, w, V) dL dw, (2.5) 
-oo JO 
roc roc ^2 

P(r) = - / / A , == =f(r,w,L)dLdw, (2.6) 



,2 



-oo JO 



y/l+w 2 +L/r 2 ' 



where p denotes the mass energy density and p is the radial pressure. 
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Before we discuss how solutions of this system can be obtained we note 
some additional equations which follow from the above. Let the quasi-local 
mass be defined by 

m(r) := 4tt / s 2 p(s) ds. (2.7) 



The ADM mass is then M = lirn r _> 00 m(r). Using m, the field equation 
together with the boundary condition (|2.1|) imply that 

e- 2 ^) = l-^, (2.8) 



p'(r) = e 2A M ( ^ + 4vrrp(r) ) . (2.9) 



and by ([23 



In particular, [i is radially increasing. Eqns. (|2,2p and f|2.3[) are the 00 and 11 
components of the Einstein field equations, and given the matter quantities 
they together with the boundary conditions suffice to determine the metric. 
But the 33 component of the Einstein equations is also non-trivial and reads 

e" 2A (// + (// + l/r)(// - A')) = Sttpt, (2-10) 

where the tangential pressure px is given by 

PT{r) = ~- / . == J (r,w,L)dLdw. (2.11) 

Eqn. (|2.10p follows from f)2.2p — (|2.6|) : the also non-trivial 44 component is a 
multiple of the former. Next we note that the Tolman-Oppenheimer-Volkoff 
equation 

P' = -(pt-p)-v'(p + p) (2.12) 
r 

holds for any sufficiently regular static solution of the spherically symmetric 
Einstein- Vlasov equation, as can be seen by a simple computation, using 
(|2.4p to express p'. Finally, if we add (j 2 . 2 j) and (|2.3p we obtain the identity 

X' + fjf = 4vrre 2A (/9 + p). (2.13) 

Let us now briefly discuss how one can establish the existence of static 
solutions of the spherically symmetric Einstein- Vlasov system; for more de- 
tails we refer to [121 13 El [T5] . Let 

E := e^ r) v/1 + w 2 + L/r 2 
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denote the local or particle energy which like L is conserved along charac- 
teristics of the static Vlasov equation (|2.4p . The ansatz 

f(r,w,L) = <i>(E,L) (2.14) 

satisfies the Vlasov equation, and the macroscopic matter quantities become 
functionals of the metric coefficient [i. Substituting these into Eqn. (|2.9|) 
and using (|2.8j) the above system reduces to a single first order equation for 
\x. Analyzing the latter constitutes an efficient way to prove the existence 
of static solutions with finite ADM mass and finite extension. It should be 
pointed out that spherically symmetric static solutions which do not globally 
have the form (12. 14ft exist, cf. [T7j. This contrasts the Newtonian case where 
all spherically symmetric static solutions have the form (|2. 14[) , a fact known 
as Jeans' Theorem, cf. [5]. As a matter of fact, below we obtain steady 
states which are good candidates for solutions which are not globally of the 
form (j2Jj]) . 

In |15l Thm. 2.1] it has been shown that in order to obtain a steady state 
of finite ADM mass by the above approach there must exist a cut-off energy 
Eq such that &(E, L) = for E > Eq and L > 0. If $ and in particular such 
a cut-off energy is prescribed we obtain a steady state solution by prescribing 
/i(0) and solving — numerically or analytically — the resulting equation (|2.9p 
radially outward. However, the resulting solution will in general not satisfy 
the boundary condition fi(oo) = 0, but it will have some finite limit n(oo). 
One can then shift both the cut-off energy and the solution by this limit to 
obtain a solution which satisfies the boundary condition /u(oo) = 0. This 
does not affect m and hence also not A, where it should be noted that the 
boundary conditions for the latter both at the centre and at infinity follow 
from (|2,8p . provided that the solution has finite ADM mass and that p is 
not too singular at the centre. 

As long as one is interested in only a single steady state at a time this 
way of handling the boundary condition for \x at infinity is acceptable, but 
it is awkward when one studies families of such steady states, since Eq and 
fi(0) cannot be treated as two free parameters if one insists on the boundary 
condition at infinity. Hence we found it convenient to use an ansatz where 
the resulting equation for [i can be rewritten in terms of the function y(r) = 
gM r ) j _ The following is sufficiently general for our purposes: 

f(r, w, L) = <f>(E/E )(L - L ) l + , (2.15) 

where I > —1/2, Eq > 0, Lq > 0, x + := max{x,0} denotes the positive part, 
and (f> :]0, oo[— » [0, oof is measurable, (f>(rj) = for rj > 1, and there exist 
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constants k > — 1 and C > such that 

m<C{l-r,)\ V G]0,l[. 

With the ansatz (pTL5]) . 

p(r) = qr 2 '(l + L /r 2 ) l+2 [g l+ - i/2 + ft+1/2 ) (^1 + Lo/r*y(r)), (2.16) 



P(r) = ^r 2 '(l + L /r 2 ) /+2 ft + 3/ 2 (Vl + L /r 2 y(r)), (2.17) 



pr(r) = (i + l)p(r) + |L r 2 ^ 2 (l + L /r 2 ) l+1 g l+1/2 (^/l + ^o/r 2 2/(r)), 
where for j > —1, Z > —1, and u G]0, 1], 

9j (u):= (j){wn){ri 2 - l) j dri = u- {2j+1) (f>(rj)(r) 2 - u 2 ) j drj, (2.18) 
^•(m) := for u > 1, and 

q := 2vr / s l (l - S y 1/2 ds. 
Jo 

The equation to be solved for y then becomes 

y'= r 4^(^ + ^ P (r,y)), (2.19) 

r 

where the dependence of m and p on y arises via the above formulas (12.161) 
and (gUP - 

In this way the cut-off energy Eq disappears as a free parameter of the 
problem, and for cj), I, and Lq fixed we obtain a one-parameter family of 
solutions parameterized by y(0) > 0. Once y is determined on the support 
of the solution its limit y(oo) at infinity can be computed, and if we define 
Eq := l/y(oo) and e^ r ^ := Eoy(r) we have a steady state with the proper 
boundary condition at infinity. 

Notice that if I = = Lq, i.e., the microscopic equation of state (|2.15|) 
does not depend on angular momentum, then = p, i.e., tangential and 
radial pressure are equal, which is why such steady states of the Einstein- 
Vlasov system are called isotropic. 

If j > -1 and k + j + 1 > then gj € C(]0, oo[). If j > and k+j>0 
then gj 6 ^QO, oo[) with 

g'j(u) = -1 ((2j + l) 9 j(u) + 2jgj_i(u)) , (2.20) 
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cf. [15\ Lemma 2.2]. With this regularity of the matter terms as functions 
of y it is easy to see that Eqn. (|2,19p has for any choice of y(0) > a unique 
solution. The non-trivial question is which choices for <f) lead to steady states 
of finite ADM mass and finite extension. In |15j Thm. 3.1] it has been shown 
that this is the case if Lq = and 

0(77) = c(l - r/) k + 0((1 - r]) k+s ) as 7] 1- (2.21) 

where c > 0, 5 > 0, and k, I € M are such that 

11 3 

fc > -1, Z > --, fc + Z + - > 0, k < I + -. 

In this case the resulting steady state is non-trivial iff y(0) G]0, 1[. 

To conclude this section we briefly discuss how we construct the steady 
states numerically. If an ansatz of the form (|2,15p is given, the value of the 
integral gj(u) and hence of p and p as functions of r and y{r) can easily be 
computed. In practice we preferred to choose ansatz functions of the form 

f(r, w, L) = (1 - E/E ) k + (L - L )+, (2.22) 

with k > 0, I > -1/2, fe < 3/ + 7/2, £ > 0, L > 0, and k, I such that the 
relevant integrals gj can be computed explicitly by hand. In [13] it has been 
shown that such an ansatz leads to finite ADM mass and compact support as 
well. In the Newtonian case with I = Lq = 0, this ansatz gives steady states 
with a polytropic equation of state. With p and p then given as functions 
of r and y(r) we use a simple Euler-type or a leap-frog scheme to solve 
the equation (I2.19|) . starting with some prescribed value y(0) g]0, 1[ at the 
centre and moving with a fixed step size Ar radially outward. If we are in a 
situation where the rigorous results cited above guarantee that the solution 
has a finite extension, the expression \/l + Lo/r 2 y(r) exceeds the threshold 
1 for sufficiently large values of r. Once this happens the computation can 
be stopped, and the spacetime can be extended by an exterior Schwarzschild 
solution of the appropriate ADM mass. 

It should be noted that when Lq > the expression + LQ/r 2 y(r) will 
always exceed 1 for sufficiently small values of r, and there will be no matter 
in the region 

r< \£S^ =:fio - <2 - 23) 

In this case y is constant in the inner vacuum region [0, Rq], and we start our 
numerical computation at the radius Rq. The choice Lq > also leads to 
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another complication, since in that case the support of the solution may in 
general consist of several concentric shells, a feature which is not captured if 
the computation is stopped the first time that y/l + Lo/r 2 y(r) again exceeds 
1. 

3 Characterization of steady states 

In this section we are interested in the possible shapes of the energy density 
p and its support. In particular we will numerically construct multi-peaked 
steady states which seem to have no analogue in the Newtonian situation. 

We start with some simple analytic observations. Without explicitly 
mentioning it we only consider non-trivial steady states of finite ADM mass 
and compact support. In particular, m{r) > for some area radius r, and 
hence m > on an interval of the form ]i?o,oo[ where Rq > is given by 
(j!T23|) . Notice that if Lq > then at r = R by definition y^/l + L /r^ = 1, 
and t-yy/i + Lo/ja = _ y (i + L /r 2 )~ 1/2 £ /r 3 < so that y y/1 + L /r 2 < 
1 on some interval of the form ]Ro, Rq + 8[. 

Observation 3.1. The function y is constant on [0, Rq] and strictly increasing 
on [Rq,oo[. This follows immediately from (12.191) . 

Observation 3.2. If the solution is isotropic, i.e., I = = Lq, then p is 
strictly decreasing and its support is an interval of the form [0, R\]. 

Observation 3.3. If Lq = 0, then the support of the solution is still an interval 
of the form [0, R\]. This follows from Observation (3.1), Eqn. (|2.16p . and 
the strict monotonicity of the functions gj on their support, cf. (12.201) . Note 
however that in the anisotropic case I ^ the energy density p is in general 
no longer monotone. 

Observation 3.4- If Lq > 0, then we have vacuum in the interval [0, Rq], 
m(r) > for r > Rq, and the support is contained in some interval [Rq, R\]. 
As we shall see below the support will in general no longer be a single 
interval. 

Steady states which are supported in [Rq, R\] with Rq > we call shells, 
and states with Rq = we call non-shells. In Section 3.3 below we give 
a quite general characterization of both shells and non-shells for different 
values of the parameters k, I, Lq and y(0). The main focus in that section is 
on cases which have rather small values on y(0) since this is a necessity for 
obtaining several peaks. In the following sections we have chosen to study 
shells and non-shells separately for values of y(0) which give no more than 
two peaks. One reason for this is that it is then easier to visualize the large 
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difference in magnitude of the first peak compared to the following one for 
these values of y(0); in Section 3.3 we really use the quantity 47rr 2 p rather 
than p itself and moreover we use a logarithmic scale in order to get more 
informative pictures for the multi-peaks. The second reason is that we want 
to emphasize the important feature that shells can have separating vacuum 
regions whereas non-shells cannot, and we find that this point is made more 
clear by first splitting the presentation into shells and non-shells. 

3.1 Shells 

Since the main purpose of this and the next section is to emphasize a couple 
of features of the shape of the energy density which are present for any 
choice k > and / > when Lq > 0, we fix k = 0, / = 1.5, and Lq = 0.2. 
A more general analysis where the influence of all parameters is taken into 
account is carried out in Section 3.3 below. 

Let us for these values of k, I, and Lq compute steady states with three 
different values of y(0), namely y(0) = 0.4, 0.22, 0.12. We find that these 

k=0, 1=1.5, L =0.2, y(0)=0.4 



1.8 - 
1.6 - 




r 



Figure 1: A pure shell 

values give rise to three distinct states. When y(0) = 0.4 we get a pure 
shell, i.e., the energy density increases, reaches a maximum value and then 
decreases to zero and remains zero for all larger values on r. In the second 
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10r 
9 - 



k=0, 1=1.5, L =0.2, y(0)=0.22 



7 - 

6 - 

o. 5 - 

4 - 

3 - 

2 - 

1 - 

- 



0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 2: A shell with a tail 



k=0, 1=1.5, L =0.2, y(0)=0.12 



55 r 

50 - 
45 - 
40 - 
35 - 
30 - 
25 - 
20 - 
15 - 
10 - 
5 - 
- 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 3: A vacuum region before the tail 
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situation, when y(0) = 0.22, the maximum value of p is followed by a local 
minimum, strictly greater than zero (this is hard to see in the figure but 
is easy to check in the data file), and then a tail with a small amplitude, 
relative to the maximum value, but with a considerable extension. In the 
final case, where y(0) = 0.12, there is a vacuum region between the first part, 
which is a pure shell, and the tail. In this case a Schwarzschild solution can 
be joined at the first point where the energy density vanishes which results 
in a different steady state, a pure shell. 

We point out that in [2 J it is indeed proved that the energy density van- 
ishes close to Ro if y(0) is small which corresponds to the third case above. 
Furthermore, in Section 2 we mentioned that Jeans' Theorem does not hold 
for the Einstein- Vlasov system, cf. [T7], and the distribution function / of 
the steady state obtained by joining a Schwarzschild solution as described 
above for y(0) = 0.12 is not globally given as one function of E and L as in 
the ansatz f)2.14j) . 

3.2 Non-shells 

In the case when Lq = 0, the energy density can be strictly positive or vanish 
at r = (depending on I) but it is always strictly positive sufficiently close 
to r = 0. Hence, the support of the matter is an interval [0, R±] with R± > 0, 
and we call such states non-shells. From Observation 3.3 it follows that no 
vacuum region can separate two parts of a non-shell. Except for this fact, 
the basic features of the shape of the energy density are preserved for non- 
shells. In Figures 4 and 5 the energy densities of two non-shells are shown, 
corresponding to k = 0, I = 1.5, y(0) = 0.4 and k = 0, I = 5.5, y(0) = 0.1. 
We see that these have very similar features as the shells in Figures 1 and 
2. Again we note (cf. Figure 5) the large difference in amplitude of the first 
peak compared to the second one. In view of Observation 3.3, which implies 
that p > in [0, Ri], it is interesting to note that non-shells nevertheless 
mimic the behaviour of shells close to r = 0, cf. Figure 5, in the sense that 
the energy density almost vanishes on an interval [0, -Ro]) f° r some small 
R > 0. 

Below we will see that the number of peaks is in general not restricted to 
two as in Figures 2, 3, and 5. However, numerically we found that already 
these solutions with two peaks seem to have no analogue for the Vlasov- 
Poisson system. 
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k=0, 1=1.5, L =0, y(0)=0.4 




1.5 



Figure 4: A non-shell without a tail 



k=0, 1=5.5, L =0, y(0)=0.1 



30 - 



25 - 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 5: A non-shell with a tail 
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3.3 Multi-peaked steady states 

In this section we construct solutions where the energy density has many 
peaks, and we give numerical evidence that there is no limit to the number of 
the peaks. Instead of plotting the energy density p itself we plot p = Airr 2 p 
(or the logarithm of this quantity) since the pictures become clearer due to 
the fact that the first peak which in amplitude is completely dominating 
gets reduced due to the factor r 2 . In general the amplitude of the peaks 
get smaller as r increases and thus the quantity 4irr 2 p shows less difference 
between the peaks and the resolution in the pictures is improved, but we 
point out that the important features of p are similar to those of p. 

We have first chosen to study the non-shells, i.e., Lq = 0. In Figure 6 the 
parameter k = is fixed and the dependence upon / and y(0) is depicted. 
We see that the parameter y(0) has an important effect on the number of the 
peaks, eg. in the last row we get from two to four peaks by decreasing y(0) 
from 0.2 to 0.05. Increasing the parameter I may also increase the number 
of peaks as seen in coloumn 2 and 3. Also the peaks become more narrow 
as / is increased. We point out that more peaks are obtained by decreasing 
y(0) further, cf. Figure 11 which shows a shell with y(0) = 0.01. 

In Figure 7 we have fixed I = 12 and y(0) = 0.05 which corresponds to 
the final case in Figure 6, and we investigate the influence of k. We find 
that by increasing k the decay rate of p is affected for large r and the peaks 
in general become wider as k is increased. 

Let us now turn to the shells. We first fix k = 1, I = 1, and Lq = 3 
and vary y(0). The result is shown in Figure 8. It is instructive to compare 
this result with the upper row in Figure 6. We see that the influence of 
a non-vanishing Lq increases the amplitude of p and also more peaks are 
obtained. The results obtained by varying Lq are then shown in Figure 9. 
First of all we notice that the number of peaks is affected, and we also note 
that the peaks get more narrow as Lq increases. In particular we get in the 
case when Lq = 1.5 one vacuum region which separates the first two peaks 
and in the case when Lq = 10, there are two vacuum regions. 

The dependence on / is shown in Figure 10 and a similar behaviour as in 
Figure 6 is found. Increasing I can give rise to new peaks and quite generally 
it makes the peaks more narrow. 

Finally we show in Figure 11 a case where fc = 1, Z = 10, Lq = 3, and 
y(0) = 0.01, which has ten peaks, and we note that there are three vacuum 
regions separating the first four peaks. 
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k=0, 1=2, y(0)=0.2 



k=0, 1=2, y(0)=0.1 



k=0, l=2, y(0)=0.05 




k=0, l=7, y(0)=0.2 



1 2 

k=0, l=7, y(0)=0.1 



1 2 

k=0, l=7, y(0)=0.05 




1 2 

k=0, 1=12, y(0)=0.2 



1 2 

k=0, 1=12, y(0)=0.1 
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k=0, 1=12, y(0)=0.05 
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Figure 6: Multi-peaks of non-shells, Lq = 
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Figure 7: Multi-peaks of non-shells, Lq = 
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k=1, 1=1, L =3, y(0)=0.35 
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Figure 8: Multi-peaks of shells 
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k=1, 1=1, L Q =0.2, y(0)=0.15 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 
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k=1 , 1=1 , L Q =1 .5, y(0)=0.15 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 



k=1, 1=1, l_ =10, y(0)=0.15 
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Figure 9: Multi-peaks of shells 
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Figure 10: Multi-peaks of shells 
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4 The Buchdahl inequality 

Let us define 

„ 2m(r) 

r:=sup — . (4.1 

r>o r 

In this section we investigate if T has an upper bound less than one for 
static solutions of the spherically symmetric Einstein- Vlasov system. We 
start with the following analytic result the proof of which is modeled on the 
original proof by Buchdahl. 

Theorem 1 The Buchdahl inequality 

r< s - 

9 

holds for spherically symmetric steady states of the Einstein- Vlasov system 
which satisfy the condition that 

p > pt and p is non-increasing. (4.2) 

These assumptions are satisfied if I = = Lq in the ansatz $2.15]) . i.e., the 
steady state is isotropic, but also if Lq =0 and I < in which case the 
steady state is anisotropic. 

Proof. First of all we notice that using (|2.13p the field equation (|2.10p can 
be rewritten in the form 

(r V~V) = 47rrV +A ( / o +p + 2p T ). 

This equation is essentially the slicing condition for maximal areal coor- 
dinates which in the static case coincide with Schwarzschild coordinates, 
cf. [I]. Using (12. 9p the above equation can be rewritten in the form 

4- ( Vy) = ^ L +2PT -1™ 3 - 2 P ) . (4.3) 

ar \r J r \ 4tt r J / 

In the sequel we abbreviate z(r) = m(r)/r 3 . Since p is non-increasing, 

. 47T ( r 2 47T f r o , 4-7T , . 

z v) = ~r P{ s ) s ds > —rp(r) / s as = — P\r). 
r 6 J r 6 J 3 

Since by assumption px < p, the function 

r 
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is non-increasing. Since 



3 47T 

z'(r) = — z(r) H p(r) < 0, 



the same is true for z. We now use 



£ = f ( r ) = 2 / se x{s) ds, r > 0, 

JO 



as a new independent variable and consider the function 

a (£) =eK»0. 

Then 

— - ' fjL [r)- 



d£ w 2re A ( r ) 

is non-increasing for r > respectively £ > 0. If we observe that Eqn. (|2.8I) 
can be rewritten as 

e A « = (1 - 2r 2 z{r))- 1 ' 2 

then for r > 0, 

a(£(r)) > a(f) - a(0) = / — (C)dC 

Jo d£ 

d£ d£ Jo 

>^(0 [ r 2s(l-2sMr))- 1/2 ds 

= e^Vfr) \, -s t~t f 1 - x/l " 2r 2 ^(r) > ) . 
P v y 2re A W z(r) V v w 7 

By (|2.9p . //(r) > e 2A rz(r), and inserting this into the previous estimate we 
find that 

e Mr) = a (£(r)) > I e MO+A(r) (1 _ e -A(0). 

This implies that e A ^ r - 1 < 3 which is equivalent to the assertion on T. 

For the ansatz (|2.15p with Lq = and I < the general assumptions 
used above are satisfied, and the proof is complete. □ 

Remark, (a) The analysis above actually yields the somewhat sharper 
estimate that 

e Mr) _ e M(0) > I e MW( e A(r) _ ^ 
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and hence 

r < 1 r-^-. 

( 3 _ 2e M0))2 

(b) Although our numerical investigation indicates that the estimate T < 8/9 
holds for steady states of the Einstein- Vlasov system in general, some key 
features used in the above proof fail, in particular, the right hand side of 
(|4.3p is not negative in general. To see this, let 

p(r) := r~ 2l p(r), p(r) := r~ 2l p(r). 

Then if I > and Lq = 0, 

r~ 21 \p + 2p T 



as r —* 0, and the right hand side of (|4.3|) is positive for r close to 0. 

(c) In [9j bounds on V are considered under various assumptions. In par- 
ticular the authors obtain the bound 8/9 under the assumptions (|4.2|) . We 
have chosen to include the above proof anyway for the sake of completeness 
and in order to stress the fact that for Vlasov matter these assumptions are 
satisfied for the ansatz (|2. 15j) with I < and Lq = so that steady states 
with these properties do exist. In [9] the authors also consider the general 
non-isotropic case without any monotonicity assumptions but instead they 
require that there is a constant B such that 

(p - 3m/r 3 ) + 2(p T - p) < 
p + m/r 3 ~ 

Under this assumption they obtain a bound on T which for large B can be 
written in the form 

T < 1 - 2/(2 + B) 2 . 

As B — > oo this bound degenerates to the estimate T < 1. This is interesting 
in view of the numerical results presented below which indicate that T is 
always less than 8/9, while the energy density along a family of steady states 
which makes T approach 8/9 gets more and more peaked so that B — > oo 
(this is easily seen by considering an approximation of a Dirac measure). 
In the case of shells such a family has indeed been constructed for Vlasov 
matter and it has been shown that T — > 8/9 in that case, cf. [2]. We also 
refer to [3] where any matter model which satisfies p+2px < £lp is considered 



4-7T r 3 



2lp 



m 



47T 7"3+2Z 



p(0) + 2lp(0) 



3 + 2/ 



p(0) > 
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(f2 = 1 for Vlasov matter) and where it is shown that a family of shells for 
which B can become arbitrary large satisfies Y < (2f2 + l) 2 — l)/(2f2) 2 so 
that r is strictly bounded away from one. Thus, as is shown in (b) above 
a proof of our general conjecture below must rely on mechanisms which are 
rather different from the ones in the proof of Theorem [T] and also different 
from [9]. 

Our next aim is to give numerical support to the conjecture that T < 8/9 
holds for all spherically symmetric static solutions of the Einstein- Vlasov 
system. Our numerical study is restricted to distribution functions of the 
form (]2.22p . Since the qualitative behaviour is quite different in the case of 
shells and non-shells we will present the results for these cases separately. 

4.1 Shells 

In Figure 12 the values of T are shown for a family of steady states param- 
eterized by y(0) € [0.03,0.6]. Here k = 0, I = 3/2, and L = 1. We see 
that r stays below 8/9. In Figure 13 we have changed the range of y(0), 
namely y(0) £ [e -14 ,e -4 ], in order to see that the family does approach 
8/9 as y(0) — > 0, which we know is the case by the result in [2]. This also 
provides a fairly tough test for our numerical scheme since the functions gj 
become very large if y gets very close to 0. 

We point out that the corresponding values of 2M/R%, where R± is the 
outer radius of the support would (as long as several peaks are present) be 
much less than T. This is clear in view of Section 3 where it was seen that 
the first peak is strongly dominating compared to the other peaks and most 
of the matter is captured in the first peak. Hence the value of r such that 
r = 2m(r)/r is within that peak. 

Let i?n be the outer radius of the first peak and consider the ratio 
Rn/Ro as y(0) -> 0. We find in Figure 14 that R u /R -» 1 as y(0) -> 0. 
This fact is indeed proved in [2] and is crucial in order to show that T —* 8/9 
for any family of shells for which y(0) — > 0. 

In contrast to the non-shell case the behaviour of V for shells is qualita- 
tively very similar for different values of k and I and we have therefore only 
included one case here. However we point out that all static shells that we 
have considered have satisfied F < 8/9. 

4.2 Non-shells 

In Figure [15] we depict T as a function of y(0) for three different families of 
non-shells, where I = 1,5,32, and k = in all cases. The results indicate 
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Figure 12: T versus y(0) for shells 



k=0, 1=1.5, L =1 
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Figure 13: T versus y(0) for shells 
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Figure 14: The ratio Rh/Rq 

that 

T < k(1) < 8/9, (4.4) 

where k is an increasing function of /. It is interesting to note the strong 
dependence of T on I which shows that the degree of anisotropy is crucial. 
In particular we see that the (non-rigorous) results by Bondi in [6] can 
be violated by anisotropic steady states. Indeed, Bondi considers isotropic 
matter models for which p > 0, p > p, or p > 3p, and gets 0.97, 0.86 
and 0.70 respectively as upper bounds of 2M/R\. Since for Vlasov matter 
p > p, we see that the second bound is violated when I is large, i.e., the 
degree of non-isotropy is large, cf. the third case in Figure [T5l For shells 
this bound is violated also when I is small as we saw in Figure 12. The 
degree of anisotropy in this case is large for any I > 0, since p/p — » as 
y(0) -0, cf 0. 

By taking I larger we get numerical evidence that also non-shells can have 
T arbitrary close to 8/9, i.e., k(1) — * 8/9, as I — > oo, but still k(1) < 8/9 for 
any I. In view of this claim and the claim made in the previous paragraph 
for shells, we find that the most important insight from our simulations in 
connection to Buchdahl's inequality is the following claim: 
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k=0, 1=1, L =0 



k=0, l=5, L =0 



k=0, l=32, L =0 




Figure 15: T versus y(0) for non-shells 



Any static solution of the Einstein- Vlasov system satisfies 

r < |. (4.5) 

As mentioned above, the result in [2] shows that the inequality is sharp 
for shells and our numerical simulations indicate that it is also sharp for 
non-shells. It is an interesting and important problem to prove (|4.5p in full 
generality. 



5 Spirals in the (R, M) diagram 

As we discussed in Section [2] for a fixed ansatz of the form (|2.15p there 
exists a one-parameter family of corresponding non-trivial static solutions 
which are parameterized by y(0) G]0, 1[ where y = e^/E^. One can now 
ask how for example the ADM mass M and the radius of the support R 
change along such a family. More specifically, one can plot for each y(0) the 
resulting values for R and M to obtain a curve which reflects how radius 
and mass are related along such a one-parameter family of steady states. As 
a first example we choose / = = Lq and k = in (|2.22p . i-e., an isotropic 
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Figure 16: (R, M) spiral for k = 0, 1 = 0, L = 0, 0.01 < y(0) < 0.95 

case, which results in the curve shown in Figure [TBI We will discuss possible 
conclusions from such curves in more detail below, but note that there are 
possibly many steady states (or none at all) which have a given ADM mass 
and differ in radius. 

As a matter of fact, for the isotropic case one can show that such radius- 
mass curves always have a spiral form. The following theorem easily follows 
from the corresponding result established by Makino for static, spherically 
symmetric solutions of the Einstein-Euler system under suitable assumptions 
on the macroscopic equation of state p = p(p), cf. Thm. 1]. For this 
theorem we parameterize our one-parameter family of steady states via the 
inverse of the central pressure, i.e., we define e := 3 , ; notice that 

g%/2 is one-to-one on ]0, 1] and e — ► iff y(0) — > 0. 

Theorem 2 Consider an isotropic ansatz of the form \2.15\) . i.e., 1 = = 
Lq, and assume that the conditions on eft specified in h2.21\) hold. Then as 

Here R(e) and M{e) denote the radius of the support and the ADM mass of 
the steady state corresponding to the parameter e, Rs and Ms are positive 
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constants, B is a non-singular matrix, 



J{9) 



cos 8 — sin 9 
sin 9 cos 9 



is a rotation by the angle —9, b S R 2 \{0} ; and 71, 72 are positive constants. 



Proof. The definition of m in (12, 7j) together with the Tolman-Oppenheimer- 
Volkoff equation (|2.12p in the isotropic case with (|2.9p and (|2.8p substituted 
in give the following system of equations: 

— — = 47rr^, (5.1) 
dr 

dp m + 4vrr 3 p 

^ = " (P+y) ^(l-2 m /r) - (5 - 2) 

Moreover, 

p = co (53/2(2/) +31/2(3/)) , 

where the right hand side is a strictly decreasing function of y G]0, 1] which 
converges to infinity as y — * 0. If we denote its inverse by 

a : [O.ooHO, 1], p^y, 

we find that the pressure p is related to the energy density p through an 
equation of state 

p = P(p):=^g 3/2 (a(p)). (5.3) 

The equations (|5.ip . (|5.2p . (|5.3p now provide a closed system of equations 
which describes an isotropic Einstein- Vlasov steady state purely in terms of 
macroscopic quantities and a macroscopic equation of state. This system 
also arises from the static, spherically symmetric Einstein-Euler system. It 
describes a gaseous star and is studied in [11]. Hence to prove the theorem 
above we only have to show that our equation of state (|5.3j) satisfies the 
assumptions made in [TT], and then we can invoke the corresponding result 
established there. 

Firstly, P = P{p) is a smooth function of p > 0, since p > iff y G]0, 1[ 
iff p > we have P(p) > for p > and since p^Oiffy^liffp^Oit 
follows that P(p) -> as p 0. Moreover, by (|2T20|> . 

dP = dP dp = 1 4g 3/2 (y) + 3g 1/2 (y) < 1 < i 

dp dy dy 3 4g 3/2 (y) + 5g l/2 (y) + g_ 1/2 (y) 3 
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for all < p < do, and also ^ > 0, so that the assumption (A.0) in [TT] is 
satisfied. 

The assumption (A.l) in [11], namely that for p — > 0, 

with some constant 4/3<7<2, is needed only to make sure that the 
solutions of the system (|5.ip . (|5.2p . (|5.3p have finite radius. Since we know 
this a priori by the assumption (|2.2ip . we need not check this condition. But 
using the results of [15] one finds that (A.l) holds with 7 = (k+5/2)/(k+3/2) 
which lies in the required interval ]4/3, 2[ iff —1/2 < k < 3/2, and this is 
precisely the assumption on k in (|2.2ip for / = 0. 

The final condition (A. 2) which we need to check is that the limit 



lim — 

p— >oo p 



a e]0,l[ 



exists, and that there exists some non-decreasing function uj such that 



P 

P 



a 



and 



1 Lo(x) 



dx < 00. 







For this it is useful to introduce the functions 



hj(y) : 



so that 



Then 



lim 



P 1 



lim 



i>(v)(v 2 -y 2 ) j dr,, y G [0,1], 

g j (y) = y~ {2j+1) h j (y), ye]0,i]. 
53/2(2/) 1 hy 2 (y) 

hm 
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P^oo p 3 y-M) 53/2 (y) + 9x/2 (y) 3 y-*0 h 3/2 (y) + y 2 ^i/2 (y) 3 ' 
note that ^3/2(0) > 0. Moreover, 



P 

P 



1 



#1/2(2/) 



lgl/2(l/) 

3 03/2 (y) +3i/2 (y) ~ 3y 3/2 (y) 

1 /c \ 1/2 51/2 (y) 1 



3 V37 (y 3/2 (y)) 1 / 2 {pogy^y)/Syi* 
1 /co \V2 yi/ 2 (y) 1 
3 V 3 



(^(y)) 1 / 2 ^' 



29 



If we can show that 51/2 (y) (53/2 (u))^ 1 ^ 2 is bounded then the condition (A. 2) 
holds with a = 1/3 and u(x) = Cy/x, where c > is a suitable constant. 
Now 

gl/aCg) _ h 1/2 (y) 

and the functions hi/ 2 and /i 3 / 2 are continuous and positive on [0, 1[. Hence 
it remains to show that the above fraction has a finite limit as y — ► 1. Let 
us denote e := 1 - y 2 -> as y -> 1. Then by [15, (3.18)], 

fygfo) _ R 1/2£ fc + 3 / 2 + Q( g fc+3 / 2+g ) ;Q 

as e — * 0, since + 3/2 > (/c + 5/2)/2 which is equivalent to > —1/2. The 
latter holds by the assumption on k in (|2.2ip . and the proof is complete. □ 
The theorem above is restricted to isotropic steady states, but numer- 
ically we have so far found no microscopic equation of state, isotropic or 
not, which does not lead to a radius-mass spiral, although for sufficiently 
anisotropic cases the spiral becomes rather deformed. Corresponding exam- 
ples are given in Figures [17] and [T8l 
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Figure 17: (R, M) spiral for k = 0, 1 = 10.5, L = 0, 0.01 < y(0) < 0.99 

The sharp corner in Figure [18] is a genuine feature of this curve and can 
be explained as follows. For Lq > the steady state can consist of several 
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Figure 18: (R,M) spiral for k = 1, Z = 5, L = 2.0, 0.01 < y(0) < 0.99 

shells of matter separated by vacuum. For y(0) close to 1 only a single 
shell is present, but at a certain value of y(0) a second shell appears at a 
certain radius which causes the outer radius of the steady state to increase 
discontinuously while the mass remains more or less constant; recall that 
p(r) > iff y(r)y/l + Lq/v 2 < 1, and although y(r) y/l + Lq/v 2 depends 
continuously on y(0) the region where the former condition holds can change 
discontinuously. 

We conjecture that such radius- mass spirals are a general feature of one- 
parameter families of steady states of the spherically symmetric Einstein- 
Vlasov system. This is surprising, since it is by no means true for the 
Vlasov-Poisson system. For example, the polytropic ansatz 

f(r,w,L) = (E -E) k + L l 

with k, I > —1, k + I > — 1/2, and k < 31 + 7/2 is known to lead to steady 
states with finite mass and compact support; the particle energy E is given 
in terms of the gravitational potential U (r) by 

E = U(r) + i(u> 2 + L/r 2 ). 

If we keep k and I fixed we obtain a one-parameter family of steady states 
which is conveniently parameterized by y(0) = Eq — U(0). The resulting 
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equation for y = Eq — U is 




47rc fcji r l y + 



21 k+l+3/2 



where Ckj is some positive constant. One can obtain any solution of this 
equation by properly rescaling a fixed, particular one, say, the one with 
y(0) = 1. Using this fact it is easy to show that ii k — I — 3/2 ^ there 
exists for each prescribed mass M e]0, oo[ exactly one steady state in the 
family with that mass, and the relation between radius and mass is 



If k — I — 3/2 = all steady states of the family have the same mass. In 
the relativistic case the picture is very different, even for the same type 
of ansatz. If the prescribed mass is too large or too small there does not 
seem to be a corresponding steady state in the family, while for masses in 
the intermediate range there are in general more than one steady state of 
that mass. For the mass M = Ms corresponding to the focus of the spiral, 
cf. Theorem [2j there are infinitely many steady states of that mass, which 
differ in their radii. 

Further interest in these radius-mass spirals comes from the fact that in 
astrophysical investigations conclusions about stability of the steady states 
in question are often drawn from them via the "Poincare turning point prin- 
ciple" : If one knows from somewhere that one steady state on the upper part 
of the spiral to the right of its first maximum is stable then it is claimed 
that all steady states up to the first maximum of M along the spiral are 
stable, and that towards the left of that point stability is lost. We can un- 
fortunately not claim to understand the arguments by which this principle is 
supported, in particular not in the case of an infinite dimensional dynamical 
system such as the Einstein- Vlasov system, and we are not aware of any 
rigorous mathematical result on stability for this system beyond the prelim- 
inary results in [18] . However, numerically we found in [4 J that along such a 
one-parameter family of steady states stability is lost if the so-called binding 
energy passes its maximum, and numerically it turns out that this happens 
at the same parameter value where M attains its maximum. Hence the 
prediction of the turning point principle agrees with the numerical stability 
analysis in [4]. 

Above we showed analytically that radius-mass spirals are not always 
present in the Newtonian case. An example where numerically one does 
find a radius-mass spiral is the so-called King model where 



M = u(k-J-3/2)/(fc+l+l/2) 



f(r,w,L) = {e E «- E -l) 
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This model is important in astrophysics. According to the turning-point 
principle outlined above there should be a transition from stable to unstable 
steady states along the corresponding one-parameter family, but it has been 
proven in [8] that all steady states of King type are non-linearly stable; the 
analogous result for the so-called relativistic Vlasov-Poisson system has been 
established in [10J. It is conceivable that in the non-relativistic case one must 
plot other quantities instead of radius and mass to apply the turning-point 
principle, but if one for example plots the radius versus the total energy for 
the King model again a spiral arises, which in view of [8] does not agree 
with the prediction of the turning-point principle. 

We conclude this section by emphasizing that the mathematically rigor- 
ous relation — if any — of such radius-mass spirals to stability properties of 
the steady states seems a very interesting and non-trivial problem in general 
and for the Einstein- Vlasov system in particular. An also very interesting 
problem is to prove that radius-mass spirals are present for the Einstein- 
Vlasov system in the anisotropic case as well, as our numerical experiments 
seem to indicate. 
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